Beyond the Hubbard-I Solution with a One-Pole Self-Energy at Half-Filling within 

the Moment Approach: Non-Linear Effects. 
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We have postulated a single pole for the self-energy, E(fe, to), looking for the consequences on the 
one-particle Green function, G(k, u>) in the Hubbard model. We find that G(k, ui) satisfies the first 
J"""""- ' two sum rules or moments of Nolting (Z. Physik 255, 25 (1972)) for any values of the two unknown 

k parameters of E(fc, u). In order to find these two parameters we have used the third and four sum 
rules of Nolting. G(k, uf) turns out to be identical to the one of Nolting (Z. Physik 225, 25 (1972)), 
which is beyond a Hubbard-I solution since satisfies four sum rules. With our proposal we have 
been able to obtain an expansion in powers of U for the self-energy (here to second order in U) . We 
present numerical results at half- filling for 1- the static spin susceptibility, x{T) vs T/t and 2- the 
band narrowing parameter, B(T) vs T/t. The two-pole Ansatz of Nolting for the one-particle Green 
function is equivalent to a single pole Ansatz for the self-energy which remains the fundamental 
quantity for more elaborated calculations when, for example, lifetime effects are included. 
Pacs numbers: 74.20.-Fg, 74.10.-z, 74.60.-w, 74.72.-h 
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After the discovery of the high-T c materials JlJ , the study of correlations has gained interested due to the fact that 
there is the belief || that the normal properties of these materials could be explained in the framework of the Hubbard 
model [|, since electron correlations are strong, i.e., the on-site electron-electron repulsions U are much larger than 
' the energies associated with the hybridization of atomic orbitals belonging to different atoms |Q . 

In this paper, we will adopt a given pole structure in the self-energy and see its effect on the one-particle Green 
function. This is done with the idea of fleshing out the basic elements implicit in the moment approach of Nolting 
for G(k,uj). At the same time, we monitor the self-energy which is the important quantity where lifetime effects can 
be incorporated rather easily. We find that by using this approach (one pole in E(fc,w) and the sum rules for the 
spectral functions), the one-particle Green function goes beyond the Hubbard-I solution, i.e., identically satisfies the 
q ' first two sum rules and its parameters have to be calculated self-consistently. 
Q . The model we study is the Hubbard model M 
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where cl (c^ ) are creation (annihilation) electron operators with spin a. = cl cj • . U is the local interaction, 
H the chemical potential and we work in the grand canonical ensemble. We have adopted Einstein convention for 
repeated indices, i.e., for the N s sites i, the z nearest-neighbor sites and for spin up and down [a = ±1). t^j = —t, 
for n.n. and zero otherwise. 

The one-particle Green function, G(k,u>), is expressed in terms of the self-energy, u>), as 

G(fc; W )EE 1 , (2) 

where e(k) — —2t(cos(k x a) + cos(k y a)), and eg = e(k) — fj,. 
We adopt the following Ansatz for E(fc,w): 

E(£,u,) = pC/+^- . (3) 

As E(fc,w) has dimensions of energy, the still unknown parameter a(k) has dimensions of (energy) 2 . a(k) is kind of 
a spectral weight and fig is the energy spectrum of the self-energy. We will calculate a(k) and fig. We have included 
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the Hartree shift directly in the self-energy (the first term of Eq. ([}])). In Eq. (||), p is the carrier concentration per spin 
orientation, i.e., p = n/2. We are assuming that we are in the paramagnetic phase ||. The physical reason of having 
chosen a single pole in the self-energy is due to the fact that, for intermediate to strong coupling, the self-energy gives 
rise to two energy branches in the one-particle and ultimately to the two Hubbard bands. 

By using Eq. (||) into Eq. (||), we get that the one-particle Green function has two poles. It can be written as 



CCl(fc) 



a 2 (k) 



LU — u>i(k) UJ — UJ 2 (k) 
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where 

u>i(k) = 

ax{k) = 
From Eqs. 



Wi(fc) 



a 2 (k) = 



4a(k) 
ui 2 (k) 
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uJ 2 (k) = - 



(%-&) 2 + 4a(fc) 



1/2" 



wi(fc) — us 2 (k) u 2 (k) - u>i(k) 

we inmediately see that the following sum rules or moments 
ai(fc) + a 2 (k) = 1 , ai(fc)wi(fc) + a 2 (k)uj 2 (k) = 



(5) 



Eqs. (|6|) are the first two sum rules for the spectral functions of Nolting 
use the next two sum rules of Nolting. This gives: 



are identically satisfied: 
^ = e(k)+pU . (6) 
J. In order to evaluate a(k) and fig, we 



u>l{k)ai{k) + u4(k)a 2 (k) = £§ + a(k) = a 2 (k) , 
uf(k) ai (k) + ul(k)a 2 {k) = e % + (Qj: + 2^)a(k) = a 3 (k) , 

where a 2 (fc), » 3 (fc) are given in Ref. |6|- as 

a 2 (k) = e 2 (fc) + 2pUe(k) + pU 2 
a 3 (k) = e 3 (k) + 3Ue 2 {k){2 + p)pU 2 e{k) + p(l - p)UB(k) + pU 3 

Temperature enters due to the presence of the B-term. Solving Eqs. (0) we find: 

a(k) = p(l - p)U 2 , n g = (1 -p)U + B(k) . 



(7) 
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The narrowing band parameter has to be calculated self-consistently and is fc-independent in the spherical approx- 
imation (See Ref. Q). Combining Eqs. (||,^) we find 



u Jl (k) = -[H(k) + [K(k)}^ 2 



u 2 {k) = - [H(k) [K(k)Y' 2 



where 



H(k) = e(k)U + B(k) , K{k) = (e(k) - U - B{k)) + 4pU (e(k) - B(k)) 



(10) 
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Eqs. ( Jill ) are nothing else that the solutions given by Nolting || by solving the four unknown cti(k) and LOi{k), with 
i = 1,2. We have regained Nolting's solutions in a much easier way, starting from the self-energy while Nolting does 
it from the Green function itself. The self-energy has the meaning of being an expansion in powers of U. In our case, 
due to the choice of a single pole for the self-energy, a(k) is of second order in U, since the first order is the Hartree 
shift, pU. 

As our solutions given by Eqs. ( [il]) satisfy the first four sum rules or moments, our Green function is beyond 
the Hubbard-I solution Q and as a consequence we have an improved solution 0. The drawback of the Hubbard-I 
solution, i.e., a gap for any value of the interaction, was pointed out by Laura Roth Q many years ago. This gap also 
exists, for all values of U ', in the spherical approximation of Nolting PJ^JlC}]. The i3-term gives rise to magnetism, a 
feature which is not present in the Hubbard-I solution. 

In Fig. 1 we present the spin susceptibility, x{T) vs T/t for the two-pole Ansatz for the one-particle Green function, 
or the one-pole Ansatz for the self-energy (See Eq. (||)), for U/t = 4.0 at half-filling. We see that the spin susceptibility 
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looks Curie-like at large temperatures and it has a bending towards zero at low temperatures. This type of behavior 
is similar to the results found for the attractive Hubbard model at half-filling. Let us remember that the attractive 
and repulsive Hubbard model become identical at half-filling. In Fig. 2 we show the band-narrow parameter, B(T) 
vs T/t for the same parameters of Fig. 1. We comment here that the band-narrowing factor, B(T), in general, is 

^-dependent. This dependence has been explicitly been shown recently to be important by Eskes and Oles O]. 

In summary we have proposed a self-energy of one pole which implies that the Green function is composed of two 
poles. This has been accomplished by the use of Dyson's equation (Eq. (||)). We have reproduced with minor effort 
the Nolting's solution ||- {§. The one-pole self-energy leading to Nolting's solution puts on firm grounds the results 
found in Refs. where it was argued that the two branches of the one-particle Green function were due to a single 
pole in the self-energy. With a single pole-ansatz for the self-energy, we recuperate known results in the literature with 
less analytical effort, at the same time sheding light on the structure of the self-energy itself. Our Green's function 
goes beyond the Hubbard-I solution because we have imposed the condition that the third and four moments be 
satisfied. For us, Hubbard-I means that the first two moments are satisfied. We add that besides these features we 
have justified that, for the particular pole structure for our self-energy, this one represents an expansion in powers of 
U. The first order is represented by the Hartree shift, or pU. (See Eq. (||)). With our rather simple approach we 
have been able to have a Green's function with two poles. Our approch is different from the one of the limit of infinite 
dimensions [ fl2"| where only the dynamical properties are taken into account leaving aside the study of the long range 
behavior. The moment approach is a reliable tool to study strongly correlated electronic systems, in particular, the 
Hubbard model. A recent calculation by Nolting, Jaya and Rex [|13| has applied it to the periodic Anderson model, 
where the relevant quantity of study is the self-energy. 
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Figures 



Fig. 1. The static spin susceptibily, x(T), vs T/t for the one-pole Ansatz for the self-energy. Here U/t — 4.0 and 
half-filling. The expression for x(T) has been borrowed from Ref. |?|. 

Fig. 2. The narrowing band parameter, B(T) vs T/t, for the same parameters of Fig. 1. 
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